Time reversal of Bose-Einstein condensates 
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Using Gross-Pitaevskii equation, we study the time reversibility of Bose-Einstein condensates 
(BEC) in kicked optical lattices, showing that inside the regime of quantum chaos the dynamics can 
be inverted from explosion to collapse. The accuracy of time reversal decreases with the increase 
of atom interactions inside BEC, until it is completely lost. Surprisingly, quantum chaos helps to 
restore time reversibility. These predictions can be tested with existing experimental setups. 
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In recent years, there has been remarkable progress in 
the manipulation and control of the dynamics of BEC 
in optical lattices (see e.g. the review In such sys- 
tems, the velocity spreading is very small. This allows 
to perform very precise investigations of the kicked rota- 
tor, known as a paradigm of quantum and classical chaos 
As a result, high order quantum resonances were re- 
cently observed experimentally Q , and various nontrivial 
effects in the kicked rotator dynamics were probed 0, S] • 
It should be stressed that in BEC the interactions be- 
tween atoms are of crucial importance, in contrast with 
other implementations of the kicked rotator with cold 
atoms [1,0, HQ, where such interactions are negligible. 
Recently, a method of time reversal for atomic matter 
waves has been proposed for the kicked rotator dynam- 
ics of noninteracting atoms [l(J. This method also al- 
lows to realize effective cooling of the atoms by a few 
orders of magnitude. The problem of time reversal of 
dynamical motion of atoms originates from the famous 
dispute between Boltzmann and Loschmidt on the origin 
of irreversible statistical behavior in time reversible sys- 
tems 11, 12 1 . The results obtained in [l(| showed that 
the quantum dynamics of noninteracting atoms can be 
reversed in time even if the corresponding classical dy- 
namics is practically irreversible due to dynamical chaos 



[13l Il4| . The investigation of the effect of interactions 
on time reversibility is of prime importance, since the 
original dispute between Boltzmann and Loschmidt con- 
cerned interacting atoms. In this paper, we study the 
effects of interactions between atoms in BEC on the time 
reversal accuracy. We emphasize that interactions bring 
new elements in the problem of time reversal compared to 
the case of one particle quantum dynamics (Io| , acoustic 
[IH and electromagnetic waves |l6| |. 

To describe the BEC dynamics in a kicked optical lat- 
tice we use the Gross-Pitaevskii equation (GPE) [l7j |: 
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where the first two terms on the r.h.s. correspond to 
the usual GPE and the last term represents the effect 
of the optical lattice, with Srit) being a periodic delta- 
function with period T. The interaction between atoms 



is quantified by the nonlinearity parameter g = NgiD 
where N is the number of atoms in the condensate and 
giD = — 2aotku± is the effective ID coupling constant, 
lu± being the radial frequency of the trap and ao the 3D 
scattering length. Here ip(x,t) is normalized to one. In 
the following, we choose units such that h = m = 1, so 
that the momentum of atoms is measured in recoil units; 
time t is measured in units of T. For g = Eq. |T]) re- 
duces to the usual kicked rotator model with the classical 
chaos parameter K = kT and effective Planck constant 
T (see e.g. |l0|). In this case the time reversal can be 
done in the way described in [Io| : the forward propaga- 
tion in time is done with T = At: + e while the backward 
propagation is performed using T — > T" = Art — e, with 
inversion of the sign of k at the moment of time reversal 
t r . This procedure allows to perform approximate time 
reversal (ATR) for atoms with small velocities inside the 
central recoil cell. This ATR procedure works quite ac- 
curately for noninteracting atoms (g — 0), however for 
real BEC the nonlinear interaction (g ^ in |T])) may 
significantly affect the return probability. Indeed, in the 
regime of strong interactions and BEC size smaller than 
the optical lattice period, earlier investigations [HI of 
soliton dynamics in GPE with g > and periodic bound- 
ary conditions showed that with good accuracy a soliton 
moves along chaotic classical trajectories of the Chirikov 
standard map As a result, even if one performs exact 
time reversal (ETR) , the presence of very small numerical 
errors completely destroys reversibility due to instability 
of chaos. However, here we are mainly interested in the 
regime typical of BEC experiments where the BEC size 
is larger than the lattice period. Our studies show that 
in this regime time reversal can be achieved with good 
accuracy for moderate strength of interactions. 

To investigate time reversibility for BEC, we numeri- 
cally simulated the wave function evolution through (JTJ) , 
using up to N s = 2 23 discretized momentum states with 
Ap = 2.5 x 1CP 5 , space discretization Ax = 2tt / (N s Ap), 
and with up to 4 x 10 4 integration time steps between two 
kicks. In this way we obtain the probability distribution 
in momentum space W p (t) = \ (p\ip(t)) \ 2 . The results of 
time reversal performed by ATR procedure are shown 
in Fig.l, for an initial Gaussian distribution in momen- 
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FIG. 1: (Color online) Initial (blue/black solid curve) and 
final return probability distributions W v obtained by ATR 
procedure vs momentum p for various nonlinearities g. Insets 
: probability distribution W p at t = t r = 10 (left) and final 
one at t = tf = 2t r (right) on a larger scale for g — 10. 
All probability distributions are scaled by their value at p = 
0, t = 0. Initial state is a Gaussian packet with rms a — 0.002. 
Here k = 4.5, T = 4tt + e for < t < t r and T" = 4tt - e for 
t r < t < 2t r with e = 2 and t r = 10. Black dots correspond 
to the wave function after ETR procedure. 

turn with rms a = 0.002 in recoil units (as in the ex- 
periment of [3j]). After this procedure, the returned wave 
packet at t = 2t r becomes clearly squeezed in momentum 
compared to the initial one. In absence of interactions 
(g = 0) the maximum of the final returned distribution 
Wq = W p =o(t = 2t r ) is equal to the maximum of the ini- 
tial one Wq = W p= o{t = 0) since ATR is exact for p = 
[Io| . The increase of nonlinearity parameter g leads to 
a reduction of the ratio Wq /Wq until complete destruc- 
tion of reversibility for g = 20. However the half-width 
Pl of the returned peak is only weakly affected by g. It 
should be stressed that at the moment of time reversal t r 
the wave packet is completely destroyed (Fig.l left inset), 
and nevertheless the peak is recreated at t — 2t r (Fig.l 
right inset). This process looks similar to the observed 
"Bosenova" explosion induced by the change of sign of 
the interactions in BEC 19]. In our case the sign of in- 
teractions is unchanged but the time reversal allows to 
invert the explosion which happens during < t < t r into 
a collapse for t r < t < 2t r . We also note that the ETR 
procedure (ip( x ) ~> ^*{ x ) at t = t r ) leads to an almost 
perfect time reversal of the wave packet, indicating that 
exponential instability is rather weak on the considered 
time scales. 

The dependence of the ratio Wq /Wq on the system 
parameters is shown in Fig. 2. It can be approximately 
described through 

Wl/Wl « e- rt " , T = Cg 2 ^/D q , (2) 

where C is a numerical constant and D q is the quantum 
diffusion rate which determines the localization length I 



FIG. 2: (Color online) Ratio of final to initial probabilities 
at p = as a function of g 2 \fat r / D q . Symbols mark the 
numerical data for: k = 4.5, a = 0.002, t T = 10 (0 < g < 20, 
blue squares); k = 8, g = 10, t r = 10 (0.002 < a < 0.2, 
green circles); k = 4.5, g = 10, a = 0.002 (1 < t r < 15, 
magenta triangles); g = 10, a = 0.002, t r = 10 (4 < k < 9, 
cyan reversed triangles); k — 4.5, g = 10 (g — 15 for the 
inset), t r = 10 (0.002 < a < 0.026, red pentagons); k = 4.5, 
g = 5 (g = 15 for the inset), a = 0.002 (1 < t r < 15, yellow 
diamonds); the curve shows the dependence ([2)| with C = 
0.22. Inset : same ratio but for a change of sign of g during 
the reversed evolution, the curve has C = 0.15. 

of quantum eigenstates in momentum space when 5 = 0: 
I — D q /2 with D q « k 2 /2 (see exact expression for D q 
in [nj). Surprisingly enough, the accuracy of time rever- 
sal increases with the increase of chaotic diffusion rate 
D q (see Fig. 2). Qualitatively, this is due to the faster 
spreading of the wave function in coordinate space, that 
leads to a decrease of the nonlinear term g\ijj(x,t)\ 2 in (JTJ) , 
making the dynamics closer to the linear case for which 
Wq — Wq. An estimate for T can be obtained assuming 
that the nonlinear term generates additional corrections 
k — * k + a(t) in ([T]), where a randomly varies in time and 
a ~ g\ip\ 2 ~ gcr/Vl- This relation takes into account the 
normalization condition (\ip\ 2 /a ~ 1) and assumes that 
inhomogeneities in 1 1/? (cc) | 2 are smoothed over / ~ D q lo- 
calized chaotic eigenstates. With these assumptions the 
decay rate of Wq /Wq is given by the Fermi golden rule 
r ~ | a | 2 ~ g 2 a 2 /D q . Such a consideration assumes that 
a remains constant during the dynamics while in fact it 
increases due to diffusion in momentum that probably 
leads to the smaller power of a found numerically 
These estimates qualitatively explain the surprising re- 
sult of Fig. [2] which shows that the increase of quantum 
chaos diffusion D q improves the time reversal at fixed g. 

In addition to the ATR procedure described above, it 
is possible to perform additionally the inversion of the 
sign of the interactions g at t = t r . In principle such an 
inversion of g can be realized in experiments similar to 
those of . The numerical data for this case are shown 
in the inset of Fig. 2. They can be described by the same 
formula as for the case of unchanged g, but with a smaller 
numerical constant C. The relatively small difference be- 
tween the two cases indicates that the main mechanism of 
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FIG. 3: (Color online) Loschmidt cooling of time reversed 
BEC atoms characterized by the ratio of final Tj and initial 
To temperatures (solid curve) as a function of k for an initial 
Gaussian momentum distribution with rms a — 0.002 and 
g = (solid curve), g — 0.5 (green/gray circles), and g — 10 
(red/black squares). Here t r = 10 and e = 2. 



time reversal destruction is related to transitions to other 
momentum states induced by the nonlinear interaction. 

The time reversal leads to a squeezing of the wave 
packet in momentum space near p — 0, that can be in- 
terpreted as an effective Loschmidt cooling [ll|. The ef- 
fect of interactions on this cooling process is analyzed in 
Fig. 3. The temperature Tf of the returned atoms can be 
defined as the temperature of atoms inside the momen- 
tum interval [—2pL, 2p^] at t = 2t r (see |l0|). For g = 0, 
the ratio 7/ /To of final and initial temperatures drops 
significantly with k. At small g = 0.5, the ratio remains 
essentially unchanged for all values of k considered. In 
contrast, for stronger nonlinearity (g = 10), there is no 
cooling at low values of k < 3, but for strong chaos with 
k > 3 the cooling reappears and becomes very close to 
the case g = at large k values. Thus surprisingly strong 
quantum chaos enhances the cooling of BEC. This result 
is the consequence of relation |(5J), according to which the 
return probability becomes larger and larger with the in- 
crease of the quantum diffusion rate D q « k 2 /2. 

The results presented in Figs. 1-3 show that time rever- 
sal can be performed for BEC through the ATR proce- 
dure even in presence of strong interactions. The cooling 
mechanism works in presence of these interactions, if the 
quantum chaos is sufficiently strong. 

Up to now we discussed the case of spatial width of 
BEC much larger than the optical lattice period. It is 
also interesting to analyze the opposite situation where 
the BEC size becomes smaller than this period. In this 
case we start from the soliton distribution 
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For k — this is the exact solution of Eq.(p}, which 
describes the propagation of a soliton with constant ve- 
locity po ■ At moderate values of k the shape of the 



FIG. 4: (Color online) Phase space distance 8 = yj Sx 2 + 8p 2 
between initial soliton and its time reversed image obtained 
by ETR (ip{x) -> if)*(x)) vs t r for k = 1, T = 2 and 
g — 10; symbols show numerical data averaged over 50 tra- 
jectories inside the chaotic (green circles) and the regular 
(red triangles) domains. The solid line shows the linear fit 
(In S) = 0.47 t r — 10.29. The slope is close to the Kolmogorov- 
Sinai entropy h « 0.45 at K — kT — 2. Insets : soliton at ini- 
tial time t = (blue/black) and final time t — 2t r (red/grey) 
with t r = 40 for initial conditions inside the chaotic domain 
(left) and regular domain (right, initial and final solitons are 
superimposed). 



soliton is only slightly perturbed, and its center follows 
the dynamics described by the Chirikov standard map 
[l8j |: po = po + fcsinxo ; xo — %o +PoT, where bars de- 
note the values of the soliton position and velocity after 
a kick iteration. In the chaotic regime with K > 1 the 
soliton dynamics becomes truly chaotic. Indeed, two soli- 
tons with slightly different initial velocities or positions 
diverge exponentially with time. As a result of this insta- 
bility, even the ETR procedure does not produce an exact 
return of the soliton to its initial state, due to the pres- 
ence of numerical integration errors. This is illustrated in 
Fig. 4, where the distance in phase space 5 between initial 
and returned solitons is shown as a function of the time 
reversal moment t r . If the initial position of the soliton 
is taken inside the chaotic domain, S grows exponentially 
with t r as 5 ~ exp(A£ r ). The numerical fit gives a value 
of A very close to the Kolmogorov-Sinai entropy h of the 
standard map at the corresponding value of K. For large 
values of t r the values of 5 become so large that the time 
reversed soliton is located far away from the initial one, 
and thus time reversibility is completely destroyed (Fig. 4 
left inset). In contrast, if the soliton starts in the regular 
domain, the growth of 5 remains weak during the whole 
integration time. In this regime, even for large values of 
t r the values of S remain small and the soliton returns 
very close to its initial state (Fig. 4 right inset). 

Another way to characterize the stability of nonlinear 
wave dynamics described by Eq.([T]) is to study the be- 
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FIG. 5: (Color online) Fidelity of two solitons with close ini- 
tial conditions vs time t for k — 1, T — 2 and different nonlin- 
earities. The initial states are solitons <(3j with g = 10 and ini- 
tial conditions (xo,po) = (0.5, 1.25) and (x' ,p' ) = (0.5, 1.255) 
inside the chaotic domain. 

havior of the fidelity defined as f(t) = \(ipi{t)\ip2(t))\ 2 , 
where ip\ and ip2 are two solitons with slightly differ- 
ent initial conditions, ft is important to note that for 
g = the function f(t) is independent of t, thus vari- 
ation with time appears only due to nonlinear effects. 
In a certain sense this quantity can be considered as a 
generalization of the usual fidelity [2l[ discussed for the 
Schrodinger evolution to the case of nonlinear evolution 
given by GPE. The dependence of / on time is shown in 
Fig. 5 for different values of the nonlinearity parameter 
g. For small values of g the function is almost constant 
on the considered time interval, while for large values of 
g it drops quickly to almost zero after a logarithmically 
short time scale corresponding to the separation of the 
two solitons. These two qualitatively different behav- 
iors can be understood as follows. For relatively weak 
nonlinearity, f(t) ~ exp(— Ft) with F oc g 2 (see Eq.©) 
that corresponds to the usual Fermi golden rule regime 
of the fidelity decay in linear quantum systems with per- 
turbation [21(. For stronger g, we enter the regime where 
nonlinear wave packets move like chaotic individual par- 
ticles that leads to an abrupt drop of fidelity as soon as 
the separation becomes larger than the size of wave pack- 
ets. In this regime the size of BEC is smaller than the 
lattice period and the time reversal is destroyed. In the 
opposite limit of large BEC size shown in Fig. [1] the time 
reversal can be maintained at moderate values of g. The 
transition between these two regimes is rather nontrivial 
and requires further investigations. 

The results of Figs. 4-5 show that the soliton dynamics 
described by GPE JT]) is truly chaotic. This leads to the 
destruction of time reversibility induced by exponential 
growth of small perturbations. However, the real BEC is 
a quantum object with a mass proportional to the num- 
ber of atoms N. Thus it has effective ft e ff oc 1/vN and 
since the Ehrenfest time t e for chaotic dynamics depends 
only logarithmically on h 0, l2l| this time remains rather 
short tE ~ \ hxh e s\/h ~ (\nN)/2h. Thus for conditions 



of Fig. 4 and N = 10 5 we have tE ~ 13 and on larger 
times quantum BEC should have no exponential insta- 
bility of motion. Hence, in the time reversal procedure 
the real quantum BEC remains stable and reversible con- 
trary to the BEC described by GPE. We think that the 
resolution of this paradox relies on the absence of second 
quantization in GPE that makes the soliton dynamics 
essentially classical. 

The dimensionless nonlinear parameter is g = 
—NaoX/2-Kl\, where l± = y / h/muj± and 2A is the op- 
tical lattice period. Values of g sa 0.1 have been achieved 
at radial frequency 800Hz and N w 5000 w N max = 
Zx/|ao|[22j. This value can be further increased up to 
g w 10 by increasing ui (e.g. 20 kHz (23|) and increasing 
A using a COvlaser [24j or tilted laser beams. Thus ex- 
perimental setups similar to 0, 0, @, l22| can test the fun- 
damental question of BEC time reversal discussed here. 
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